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Abstract 



Glucose homeostasis is controlled by the islets of Langerhans which are 
equipped with a-cells increasing the blood glucose level, /3-cells decreasing 
it, and 5-cells the precise role of which still needs identifying. Although 
intercellular communications between these endocrine cells have recently 
been observed, their roles in glucose homeostasis have not been clearly un- 
derstood. In this study, we construct a mathematical model for an islet con- 
sisting of two-state a-, (3-, and 5-cells, and analyze effects of known chemical 
interactions between them with emphasis on the combined effects of those 
interactions. In particular, such features as paracrine signals of neighboring 
cells and cell-to-cell variations in response to external glucose concentrations 
as well as glucose dynamics, depending on insulin and glucagon hormone, 
are considered explicitly Our model predicts three possible benefits of the 
cell-to-cell interactions: First, the asymmetric interaction between a- and 
/3-cells contributes to the dynamic stability while the perturbed glucose level 
recovers to the normal level. Second, the inhibitory interactions of 5-cells 
for glucagon and insulin secretion prevent the wasteful co-secretion of them 
at the normal glucose level. Finally, the glucose dose-responses of insulin 
secretion is modified to become more pronounced at high glucose levels due 
to the inhibition by 5-cells. It is thus concluded that the intercellular com- 
munications in islets of Langerhans should contribute to the effective control 
of glucose homeostasis. 
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1 Introduction 



Homeostasis, maintenance of the constant physiological state, is one of the 
main characteristics of life. In particular, glucose homeostasis is critical 
because glucose is the energy source for our bodies; the malfunctioning of 
this process causes several disease states including diabetes mellitus and 
brain coma. 

In order to understand glucose homeostasis, we first need to examine 
the tissue controlling the blood glucose level, the islet of Langerhans in the 
pancreas. It consists mainly of three types of endocrine cells: a-cells which 
secrete glucagon hormone increasing the glucose level, /3-cells which secrete 
insulin decreasing the glucose level, and 5-cells which secrete somatostatin, 
known to inhibit activities of a- and /3-cells. The hormone secretion of a cell 
influences the behavior of neighboring cells, and is thus tightly correlated 



with the islet structure (H opcroft et al., 1985[ |Pipeleers et al., 1982[ ). In 
rodents, an islet contains about 1,000 endocrine cells on average: /?-cells, 
occupying the most volume (70 to 80%) of an islet, populate largely in its 



core, whereas non-/3-cells are located on the mantle (Brissova et al., 2005). 

To the first approximation, a- and /3-cells should be sufficient for glucose 
control because a-cells can increase the glucose level whereas /3-cells can de- 
crease the level. The importance of this bi-hormonal mechanism for glucose 



homeostasis has been well recognized (Cherrington et al., 1976). However, 
it should be noted that endocrine cells in the islet interact with each other 
rather than act independently. For example, the electrical coupling between 
/3-cells through gap-junctions is known to enhance insulin secretion of cou- 
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pled /3-cells flJo et al, 2005] |Pipeleers et al., 1982||Sherman et al, 1988P - In 
addition, it has been recently reported that chemical interactions between 



neighboring cells through hormones (Cherrington et al., 1976; Franklin et al., 2005 



Orci fc Unger, 1975 


Ravier fc Rutter, 2005 


Samols et al, 1965 


Samols fc Harrison, 1976 



Soria et al, 2000) and neurotransmitters (Brice et al., 2002 Franklin fc Wollheim, 2004 



Gilon et al, 1991t[Mo7iyama fc Hayashi, 2003[|Rorsman et al., 1989[|Wendt et al, 2004Q , 
termed "paracrine interaction," affect glucose regulation. 

Among these intercellular communications, enhancement of insulin secre- 



tion by glucagon ( Brereton et al., 2007 Samols et al, 1965 Soria et al., 2 000) 
seems to be paradoxical because a-cells, playing the reciprocal role to /3-cells 
in glucose regulation, promote the activity of /3-cells. In contrast, insulin, se- 



creted by /3-cells, inhibits glucagon secretion of a-cells (Cherrington et al, 1976 



Franklin et al, 20051 IRavier fc Rutter, 20051 ISamols fc Harrison, 19761 ISoria et al, 2000p, 



which appears natural. Furthermore, the role of the third cell-type, 5- 
cells, is still not completely known although there have been reports that 
somatostatin hormone, secreted by 5-cells, suppresses the hormone secre- 



tion of both oi- and /3-cells (Cherrington et al., 1976 Daunt et al., 2006 



Orci fc Unger, 1975| [Soria et al, 20001 ) . 

What is then the raison d'etre of the paradoxical interactions between 
a- and /?-cells and the inhibitory action by (5-cells? Despite previous stud- 



ies as to these questions over the last thirty years (Orci fc Unger, 1975 



Pipeleers, 1987 ; Soria et al., 2000 Unger fc Orci, 1977 ) , there still lacks con- 
crete understanding of the role of these interactions in terms of glucose 
homeostasis. The primary difficulty in understanding these interactions lies 
in the complexity of the islet system which includes many interactions be- 



tween different coexisting cell- types working in different conditions. 

In this paper, we analyze the interactions between a-, (3-, and <5-cells, 
which contribute to the precise control of the glucose level, by means of a 
mathematical model incorporating experimentally known interactions be- 
tween islet cells (see above). As a result, our model predicts that the in- 
tracellular interactions modify insulin and glucagon secretion in a way to 
control the blood glucose level more efficiently. 

2 Islet model 

2.1 Activity of islet cells 

We begin with a simplified model in which cells of each type can take one 
of two states (active and silent). The state of a cell, represented by "Ising 
spin" a, is defined to be active (a = +1) when the cell secretes islet hormone; 
otherwise the state is defined as silent (a = — 1). Accordingly, the state of 
an islet consisting of a-, (3-, and <5-cells can be represented by {a a ,ap,as). 
There is a total of 2 3 possible states of the islet, among which (+1, —1, —1) 
and (—1,-1-1,-1-1) describe the islet state at low and high glucose levels, re- 
spectively. The main source of changing cell states is the blood glucose level 
G, which globally influences all cells. In addition, the paracrine interaction 
J from neighboring cells locally affects cell states. In this manner we obtain 
a simple Ising-type model, generally characterizing two-state dynamics in 
statistical physics: The glucose level G corresponds to the external mag- 
netic field and the paracrine interaction J to the local interaction between 
spins. 
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,iFrom the known cellular interactions illustrated in Fig. [IJ one may 
determine local stimuli G a , Gg, and Gs, which change the states of a-, (3-, 
and 5-cells, respectively, in the forms: 



n - r l + a P j 1 + <7 <* t 

<^a — — <J ^ a P 2 S 



Ug — GH Jp a Jgs 

Gs = mG, 



where G = G — Go measures the excess glucose level from the basal glucose 
level Go during the fasting period. The reciprocal nature of a- and /3-cells 
in the responses to glucose is manifested by the opposite signs in front of 
G in the first equation (for G a ) and the second one (for Gg) of Eq. [TJ In 
addition, the asymmetric interaction between these two cell types is also 
reflected in the second terms involving J a g and Jg a of the equations. For 
simplicity, we assume that the interaction strength Jg a from a- to /3-cells is 
the same as J a g from (3- to a-cells and given by Ji, i.e., J a g = Jg a = J\. 
The last terms involving J a s and Jgs describe the inhibition effects of 5- 
cells on a- and /3-cells, both with negative signs. Although the endogenous 
strengths of the interactions from <5-cells to a- and /3-cells are not known, 
the exogenous stimulus of somatostatin has been reported to inhibit both 



insulin and glucagon secretion to a similar degree (Cherrington et al, 1976). 
As a first approximation, it is thus assumed that both interactions have the 
same strength: J a s = Jgs = Ji- Here the interaction strengths J\ and 
J2 are expressed in terms of the relative effects to glucose stimulation, and 
therefore have the unit of mM corresponding to the hormonal stimulus J, 



5 



namely, a given amount of stimulus J by hormone is considered to produce 
the same effects on a cell as a certain amount J of glucose stimulation. In our 
simplified model, 5-cells are not influenced by neighboring a- and /3-cells but 
stimulated solely by glucose; therefore, Gs depends only on G in Eq. [TJ Like 
/3-cells, (5-cells become active, and secrete somatostatin above a threshold 
level of glucose. The glucose sensitivity of <5-cells is expected to have a value 
between zero and unity, i.e., < m < 1 because the threshold level for 



the activation of <5-cells is lower than that of /3-cells (Efendic et al., 1979 



Nadal et al., 1999 ). We thus choose the value m = 0.5 in this study; the 
overall behavior does not depend qualitatively on the value of m. 

For given local stimulus G a considering glucose stimulus G and effects 
of insulin and somatostatin, the transition rate of an a-cell from state a a 
to state —a a depends on the states of other cell types as well as its own 
state, and is denoted as w a {a a , ap, as). This transition rate should satisfy 
the detailed balance condition between two a-cell states a a and —a a at 
equilibrium: 

w a (a a ,ap,as)P{a a ) = w a (-a a ,a/3,as)P(-a a ), (2) 

where the probability P{a a ) for state a a follows the Boltzmann distribution 
exp[— G a (l + <7 Q )/2@] with respect to the quantity G a (l + cr a )/2 for the 
local stimulus G a . Namely, a-cells favor the state minimizing the quantity 
G a (l + a a )/2. Here measures the amount of uncertainty, which is in- 
evitable in biological systems. The origin may be the heterogeneous glucose 
sensitivity of cells and/or the environmental noise including thermal fluc- 
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tuations. It is obvious that G a (l + a a )/2 and correspond to the energy 
and the temperature, respectively, in statistical physics. In this study, the 
"temperature" is taken to be unity (0 = 1) in units of the "energy", which 
is biologically tantamount to the fluctuations caused by 1 mM change of 
glucose stimulation. 

The ratio between the reciprocal transition rates thus reads 



w a {-a a ,a/3,as) 



exp 



cxp 



— G a a 



q I G e a a a + -^OuGp + -j-a 5 <J a 



Ji 



J 2 



(3) 



with Gf = G + Ji/2 + J 2 /2, where Eq. Q] has been used to obtain the 
second line. There the three stimulation terms represent effective glucose 
stimulation, paracrine interaction from /3-cells, and another from <5-cells, 
respectively. Assuming that these stimuli affect independently the ct-cell 
state, we write the transition rate in the form 



1 

x 



1 + tanh 



G\ 



wr 



20 



I + tanh ( ^ ) a„<7 s 



I + tan]] ( ^ | a^rr,, 



(4) 



where r measures the characteristic time of the transition and it has been 
noted that tanh(yo") = a tanh y for a = ±1. Note that among possible tran- 



sition rates satisfying Eq.[2l we adopt the Glauber dynamics (Glauber, 1963) 
to choose the specific form of Eq. [H which exhibits the sigmoidal form ubiq- 
uitously describing response functions in biological systems. However, the 
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behavior of the system in general does not depend qualitatively on the spe- 
cific form of the transition rate satisfying Eq. [2j Similarly, we obtain the 
transition rates wp and w$ of (3- and <5-cells. 

The transition rates of three cell types can be summarized as 

w x (a a ,ap, a s ) = -!- [w x + w x a a a + w x pOp + tof a s + w x a pa a ap 

+Wp S apa s + w x a a 5 a a + w%p 5 a a apa s ] , (5) 

with x = a,/3, and 5, where the coefficients are given in Tables [1] and El The 
master equation, describing the evolution of the probability P(a a , ap, as) 
for the islet in state (a a , up, as), reads 

—P(a a ,ap,as) 

= w a (-a a , op, a s )P(-a a ,ap, a s ) + wp(a a , -ap, a s )P(a a , -op, a s ) 
+ w$(a a , ap, -a s )P(a a , ap, -a s ) - w a (a a , ap, a$)P(a a , ap, a s ) 
- wp(a a , ap, as)P(a a , ap, as) - ws(a a , ap, as)P(a a ,ap, as) (6) 

with the transition rates w a , wp, and wg in Eq. Note that Eq. [U] describes 
the net flux to state (a a , ap, as) simply given by the difference between the 
in-flux to state (a a , ap, as) from other states (— a a , ap, as), (a a , —ap, as), 
and (a a ,ap, —as) and the out-flux from state (a a , ap, as) to others. 

^From this master equation, it is straightforward to obtain the time 
evolution of the ensemble averages of the cell states and their correlations. 
For example, multiplying both sides of Eq. [6] by a a and summing over all 
configurations, we obtain the evolution equation for the average (a a ) = 



S 



cr a P(cr a ,(jp,as) of the state of a-cells: 




2(cr a tW Q: (cr a , as)) 



(7) 



and similarly, 




2(apwp(a a ,ap,as)) 



2(asw a (a a ,ap,as)). 



(8) 



0) 



Note that 1 + (a a ) gives twice the average activity of a-cells, etc. 

Among the eight equations for the probability P(a a , a/3, as) correspond- 
ing to the eight possible states of the islet, only seven are independent, due 
to the normalization condition ^2 aa a as p(^a,^/3,crs) = 1. Therefore, there 
exist four more equations in addition to the above three describing the av- 
erage of cell states. Those are evolution equations for correlations of two 
cell states and of three cell states. The equation for the correlation func- 
tion (a a ap) of the a-cell and /3-cell states can again be derived from Eq. [6l 
multiplied by a a ap and summed over all configurations: 

= -2(a a Of3W a (o a ,Of3,os)) - 2{a a a/3 wp (a a , 073, as)) . (10) 

The equations for (apas) and (asa a ) are also obtained in the same way: 

^( a l3 a s) = -2(vf3V5Wp(cr a ,ap,a s )) -2{apa 5 ws(a a ,a(3,as)) 
^((rsVa) = -2{a 5 a a w 5 {a a ,a l3 ,as)) -2{asa ol w a (a a ,a l s,as)). (11) 
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Note that correlations of two cell states represent the relative activity of the 
two cells. Accordingly, it makes a good measure of the different responses 
between two cells. Similarly, the equation for correlations of three cell states 
is given by 

-2(a a ai3asw s (a a ,ai3,as)). (12) 
Substituting the transition rates in Eq. [5] into Eqs. [7] to [T2l we finally 
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obtain equations for the states of the three cell types and their correlations: 

T J t ( a <*) = ~K - w a {a a ) - wf^ap) - wf a (a 5 ) - w%(a a ap) 
- w Zps( a P a s) ~ wf(a s a a ) + Wp S (a a apa s ) 

-w%(apa s } - w^ s (a s a a ) - w^ a {a a apa s ) 

-w 5 p{apa 5 ) - w s a (a S (J a ) - w & a p(a a ap(T5) 

-(wf a + w^Jiapas) - (w% s + w% s )(a s a a ) 
-(wf + w^)(a a apa 5 ) 

-(w 13 + w s ){a p a 5 ) - (w% + w 5 afS )(a s a a ) 
-{w P a + w S a ){cr a <jpo 5 ) 

Tj t {G5°a) = -(wj a + Wfj ~ {wj + wf )(<T a ) ~ {w af}5 + W^ S )((Tp) 

-(w s a/3 + w%p){<Tp<T S ) - (w s + w a )(a 5 a a ) 

T^i^crpcrs) = -(v%ps + w^pg + w S aj35 ) - (wfj 5 + w% 5 + w S ps )(a a ) 

-(wf + w% + w 5 5 )(a a a p ) - (wf + tfi + w s a )(apa s ) 
-{w% + w% + wffl(a s a a ) - (w a + w 13 + w/ 5 )<<7 a <7 /3 a ( 0L3) 



2.2 Glucose homeostasis 

Heretofore we have focused on the cellular interactions at a given glucose 
level. To study dynamics of glucose homeostasis, however, we should also 
take into account the change of the glucose level and incorporate another 
equation for glucose regulation into the model. Based on the fact that a- 
and /3-cells secrete glucagon and insulin, respectively, raising and reducing 
the glucose level, the equation for the glucose level G is taken to be 

dG _ 1 + (a a ) _ 1 + M . , 

G dt 2 2 ' 1 ' 

where tq is the characteristic time for the hormones to regulate the glucose 
level. It is expected that tq is larger than the characteristic time r of the 
change in cell states. Equation [14] describes the decrease or increase of the 
glucose level when a-cells or /3-cells are active {a a = 1 or ap = 1). Here, 
for simplicity, we have used the same characteristic time tq for glucagon 
and insulin to regulate glucose levels. Having different time constants turns 
out merely to shift the stationary level of blood glucose. To sum, we have a 
total of eight differential equations given by Eqs. [13] and [21 which describe 
the process of glucose homeostasis. 

3 Results 

3.1 Asymmetric interactions between a- and /?-cells 

In our model, activities of a-, (3-, and 5-cells are determined by the exter- 
nal glucose level together with feedback loops of intercellular interactions. 
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A given cell, subject to a glucose stimulus, secretes hormone which influ- 
ences the behavior of neighboring cells. In response, the neighboring cells 
reversely influence the given cell. These mutual interactions through hor- 
mones constitute the feedback loop which is widely employed for advanced 



system control in engineering (Bechhoefer, 2005) 



The interactions between a- and /3-cells are asymmetric: While glucagon 



secreted from a-cells enhances insulin secretion of /3-cells ( Brereton et al. , 2007 



Samols et al, 1965; Soria et al., 2000 ), insulin inhibits glucagon secretion (Cherrington et al, 1976 



Franklin et al, 2005; 


Ravier & Rutter, 2005 


Samols & Harrison, 1976 


Soria et al, 2000) 



The former positive interaction to the counterpart cells may seem strange, 
but it eventually contributes to the construction of a negative feedback loop 
for both cells. At low glucose levels, a-cells secrete glucagon, which enhances 
insulin secretion. In turn, insulin inhibits the glucagon secretion of a-cells. 
Therefore, their interactions as a whole tend to suppress the glucagon se- 
cretion from a-cells. Similar negative feedback operates when /3-cells are 
activated by high glucose concentration. It is noteworthy that this feedback 
works more efficiently in case that the glucose level varies. At a static glu- 
cose level, it should be difficult for the mutual interactions between a- and 
/3-cells to arise simultaneously because insulin and glucagon are secreted at 
different glucose levels. 

In general, a negative feedback favors stability of a system because it 
attenuates overaction of the system such as overshoot or undershoot. The 
negative feedbacks in an islet system contribute to the stable recovery to 
the normal glucose level Goo, when the system is externally perturbed by 
stimuli such as a glucose dose. The normal glucose level Goo, reached by 
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G (= G — Go) at stationarity, depends on the cellular interactions shown in 
Fig. [2j The asymmetric interaction J\ lowers the basal glucose level because 
a-cells activate /3-cells which secret insulin and thus reduces the glucose 
level. In addition, the inhibitory interaction J2 of <5-cells, albeit the same 
for a- and /3-cells, suppresses the activity of /3-cells more than that of a-cells 
at the normal glucose level, because the activity of /3-cells is higher than that 
of a-cells resulting from the asymmetric interaction between a- and /3-cells. 
Accordingly, the basal glucose level tends to increase as the strength J2 of 
the inhibitory interaction is increased. 

Figure [3] demonstrates the smooth recovery of the glucose level in the 
presence of cellular interactions (solid line), compared with the somewhat 
erratic recovery, once reaching low glucose levels, in the absence of the inter- 
actions (dashed line). For comparison, we also consider the behavior in the 
case of symmetric interactions between a- and /3-cells, i.e., where glucagon 
inhibits insulin secretion and vice versa, only to find even more erratic recov- 
ery (see the dotted line). Shown here is the recovery from the high glucose 
state [G = ImM (or G = G + ImM), (a a ) = -1, and (073) = 1]. The 
recovery from a low glucose state gives the same results (data not shown) 
although such erratic recovery is more pronounced for the glucose level start- 
ing from a higher value. 

To examine the stability in approaching the normal glucose level, we 
define the balance function 



6(G) S r/4- = i±<^ - i±M, (15 ) 
K ' u dt 2 2 y J 
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which describes the glucose level change during the characteristic time. Since 
the activity of cells represents their hormone secretion, b(G) appropriately 
describes the effectiveness of the glucose regulation by a- and /3-cells. If 
the characteristic time tq of glucose regulation is much larger than the 
characteristic time r of cell responses in Eq. [I3j i.e., r <C tq, the glucose 
level should be in a quasi-stationary state at time t shorter than tq. Then the 
fast dynamics of cell states in Eq. [13] saturates rapidly at a given glucose level 
and the seven variables, activities and correlations, reach their fixed points 
depending on the glucose level G. In particular (a a ) and (erg) depend on G, 
giving the balance function in Eq.[15]as a function of G, with a fixed point at 
G = Goo (see Fig. [4]). At low glucose levels (G < Goo), we have the balance 
function greater than zero (b > 0), or dG/dt > 0, thus the glucose level 
grows with time. At high glucose levels (G > Goo), the opposite behavior 
arises. The resulting flow of the balance function is illustrated by the arrows 
in Fig. S] and it is concluded that the balance function correctly describes 
glucose homeostasis. Further, the slope of b(G) near the fixed point G — Goo 
represents how smoothly the glucose level approaches the normal level: The 
slope of the balance function for the asymmetric interaction is small at 
G = Goo, which results in the smooth recovery of the normal glucose level 
shown in Fig. O This result is more evident with the interaction strength 
J\ larger and the characteristic time tq shorter. 

If a- and /3-cells would inhibit each other, how should the result change? 



As suggested already (Saunders et al, 1998), the bidirectional inhibitory in 



teractions seem to be optimal in view of that a- and /3-cells play opposite 
roles in glucose regulation. Remarkably, however, such symmetric interac- 
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tions turn out to result in dynamically unstable responses, as shown by the 
dotted line in Fig. (3} If this were the case, glucagon secreted by a-cells 
at low glucose levels would suppress /3-cells from secreting insulin. As the 
secretion of insulin decreases, so would the inhibitory effects of insulin on 
the glucagon secretion diminish. It should thus follow that glucagon secre- 
tion is not negatively controlled, implying more glucagon secretion. Such an 
apparent positive feedback loop, enhancing hormone secretion, gives rise to 
an instability in the islet system (see Fig. [3]). 

3.2 Inhibitory interactions of 5-cells 

3.2.1 Suppression of co-secretion from a- and /3-cells 

There is basal hormone secretion from a- and /3-cells even at the normal 



glucose level (Cherrington et al., 1976), where it is not necessary to change 



the blood glucose concentration with the help of glucagon or insulin. Ob- 
viously, the simultaneous secretion of glucagon and insulin at the normal 
level should be minimized because the opposite effects of the two would 
cancel out, nullifying the net effects on the glucose level. Such wasteful co- 
secretion of counteracting hormones can be prevented by <5-cells secreting 
somatostatin, which inhibits secretion of both glucagon and insulin. 

In our model, the average activity of cells is given by (l + (<r))/2. Accord- 
ingly, the average cell state (a) = ±1 means that all cells are active/silent; 
in particular (er) = corresponds to half of the cells being active. In the 
absence of the inhibitory interaction of (5-cells, Fig. EJa) shows that both 
(cT Q } and (an) take values greater than —1 even at the normal glucose level. 
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Namely, fluctuations associated with the biological uncertainty have some 
fraction of cells still active, leading to basal hormone secretion. Here the 
presence of inhibitory interactions of <5-cells lowers the basal activity of a- 
and /3-cells, as shown in Fig. [5jb), which reduces co-secretion of the coun- 
teracting hormones, glucagon and insulin. 

Figure O displays the relation between (a a ) and (<Jg), in the absence 
(J2 = OmM) and presence (J2 = 2mM) of the inhibitory interaction of 5- 
cells. The system at low or high glucose levels is described by the upper 
left or lower right parts of the curves on the ((073), {cr a )) plane, respectively. 
Namely, when the glucose concentration is low, a- and /3-cells are in high 
and in low activity, respectively ((cr a ) > and (073) < 0); this is reversed at 
high glucose concentrations. It is manifested that the inhibitory interaction 
of (5-cells reduces simultaneous activation of a- and /3-cells. Compared with 
the result for J2 = mM (dashed line) , the result for J2 = 2 mM (solid line) 
shows that the activity of (3- or a-cells is reduced substantially at high or 
low glucose levels. In particular a-cells remain almost silent ((<r a ) ~ —1) at 
high glucose levels. Note, however, that those endocrine cells are not totally 
silent at given glucose levels and still exhibit residual activity, which results 
from fluctuations in the glucose responses of the cells. Interestingly, it was 
suggested that such basal hormone secretion also plays an effective role: 
The minimal basal secretion of glucagon compensates the glucose uptake in 
the liver while basal secretion of insulin inhibits over-secretion of the basal 



glucagon (Cherrington et al, 1976). 
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3.2.2 Enhancement of glucose dose-responses of /3-cells 

Another consequence of the inhibitory interaction of 5-cells is the shift of 
glucose dose-responses for insulin secretion to the right direction. This is 
associated with the increased control of /3-cells by 5-cells at high glucose lev- 
els. Figure [7] indeed shows that the shift leads to more conspicuous glucose 
responses of /3-cells at high glucose levels. In general /3-cells are coupled 
with each other through gap-junction channels, which help the cells syn- 



chronize their behaviors (Sherman & Rinzel, 1991). A /5-cell cluster thus 



tends to produce all-or-none glucose responses ( Soria et al, 2000 ). In the 
real islet, on the other hand, <5-cells, with their inhibitory interactions de- 
pending on the glucose level, can modify the glucose dose-response of /3-cells. 
Accordingly, insulin response can be more pronounced at high glucose levels 
(G > 0). 

It is observed that some primitive animals have only (5- and 5-cells in 
their islets, unlike the mammals whose islets contain a-cells as well as (3- 



and 5-cells QFalkmer, 1985[ >. This difference could perhaps be attributed to 
an evolutionary adaptation. At early evolutionary stages, the islet might 
be a passive system: Without a-cells directly increasing the glucose level, 
the glucose level should increase passively as a result of the decrease in 
insulin secretion. Still, the precise glucose dose-responses at high glucose 
levels could be possible with <5-cells. At later stages, equipped with a- 
cells, the islet became an active system with regard to glucose regulation. 
It is of interest that this evolutionary change is correlated with the fact 
that (3- and 5-cells are closer to each other than a-cells in the develop- 
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ment of a stem cell (Kemp et al., 2003). In addition, 0- and <5-cells have 
functional similarities of using ATP-dependent K + channels in glucose re- 



sponses (Quesada et al, 2006 Quesada et al, 1999) 



4 Discussions 

The islet of Langerhans is a precise system that controls the glucose level 
through the use of three main types of endocrine cells. Here it is of in- 
terest to investigate whether the existing interactions between those cells 
are beneficial for glucose homeostasis. There are some evidence for the 
critical role of the interactions, which may not obviously be addressed by 
probing a- and /3-cells separately. The molecular mechanism of how a-cells 
regulate glucagon secretion at variable glucose levels is still not clearly un- 



derstood ( |Gromada et al, 2007[ ). Several works attempted to explain this 
by means of the interactions between a- and /3-cells: At high glucose lev- 
els, glucagon secretion is inhibited by insulin, GABA, or Zn 2+ secreted from 



/3-cells ( |Gromada et al, 2007[ [Ishihara et al., 2003[ ). There is also a hypoth- 
esis that glucose has direct effects on a-cells through endoplasmic reticulum 



Ca + storage (Vieira et al, 2006). Another evidence for the role of cellular 



interactions in glucose homeostasis comes from hyperglucagonomia, which 
occurs in diabetics at abnormally high glucose levels. It appears paradoxi- 
cal that the glucagon levels of such patients are high even though the blood 



glucose levels are high enough to make a-cells silent ( |Gromada et al., 2007[ ) . 
This puzzling result can be explained on the basis of cellular interactions in 



an islet (Franklin et al, 2005; Rorsman et al., 1989; Takahashi et al., 2006). 
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Note that there is also another explanation of this phenomenon in terms of 



the peculiar glucose dose-responses of (rat) a-cells ( Kemp et al, 2 003). 

In contrast, there also exist a few reports that some cellular interactions 
may not exist and are not necessary for glucose homeostasis: It has been 
proposed that the microcirculation from (3- to peripheral a- and (5-cells pro- 



hibits the paracrine action from non-/3 to /?-cells (Wayland, 1997). In addi 



tion, it has recently been reported that islet transplantation is successful in 



recovering from hyperglycaemia with only /3-cell clusters (King et al., 2007). 

Nevertheless, the existence of the receptors of signalling molecules such 
as insulin, glucagon, somatostatin, glutamate, and GABA, which are ex- 
pressed in pancreatic endocrine cells, apparently implicates their physi- 



ological roles in the fine control of glucose levels (Gromada et al., 2007 
Strowski & Blake, 2008). A better understanding of this tissue, therefore, 
will contribute to more advanced medical treatment of diabetes than the 
current one relying mostly on insulin. For example, it is conceivable to 
use other hormones such as glucagon and somatostatin for more active and 
precise glucose control. 

A variety of complicated interactions in an islet makes it difficult to 
recognize their roles, and existing experiments as to those interactions have 
focused mostly on static responses of the endocrine cells. However, it is likely 
that the cellular interactions actually contribute to dynamical responses to 
glucose. In this study, therefore, to understand the role of intercellular in- 
teractions between a-, (3-, and 5-cells, we have proposed an islet model and 
investigated the effects of integrated intercellular communications between 
those cells in glucose homeostasis. Our mathematical model can systemati- 
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cally include all the cellular interactions and identify their effects on static 
and/or dynamic responses to external glucose changes. It also takes individ- 
ual heterogeneities into consideration, e.g., in glucose sensitivity; the basal 
hormone secretion at the normal glucose level reflects that some cells can be 
active to secrete hormones even though most of the cells are silent at that 
glucose level. The small variations in glucose responses among homologous 
cells may contribute crucially to the cellular interactions between heterol- 
ogous cells, which are actually activated in quite different glucose concen- 
trations, because the heterogeneous responses of homologous cells can lead 
to an overlap in the activation between the heterologous cells. ^From this 
model, it has been revealed that the interactions give more stable, efficient, 
and accurate control of glucose: First, asymmetric interactions between a- 
and /3-cells contribute to the dynamic stability when the glucose level, per- 
turbed from the normal level, recovers to the latter. Second, the interactions 
of somatostatin for glucagon and insulin secretion prevent their wasteful co- 
secretion at the normal glucose level. In addition, at high glucose levels, the 
inhibition by <5-cells modifies glucose dose-responses of insulin secretion. For 
a more realistic and accurate understanding, it would be necessary to know 
the physiological values of the model parameters. In particular, the relative 
effects of the direct glucose stimulus G and paracrine interactions J on the 
states of endocrine cells should be identified. 

Here it is proposed that these predictions can be verified in experiment. 
As for the role of 5-cells, our results may be confirmed through the use 
of cell clusters of different compositions of cell-types, for which the culture 



method was used in the existing study (Pipeleers et al, 1982). Another 
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prediction related with the asymmetric interactions between a- and /3-cells 
needs to be verified in vivo experiment on transgenic mice because the effects 
should arise in the dynamics of whole-body glucose regulation. Note that 
the specific cellular interaction may be blocked selectively in knockout mice 



lacking specific hormone receptors in an endocrine cell (Diao et al., 2005 



Sorensen et al., 2006). 



Beyond the interactions between endocrine cells analyzed in this study, 
there exist reports that <5-cells are also influenced by a- and /3-cells ( Unger &; Orci, 1977 ) 
and these paracrine interactions should be considered with the microcircula- 



tion of hormones in an islet as described above (Wayland, 1997). It has also 
been reported that there exist autocrine interactions via which a cell is af- 
fected by its own hormone secretion ( Aspinwall et al., 1999 Cabrera et al., 2008 ). 
Furthermore, input from exocrine cells (Bertelli et al., 2001 Bishop & Polak, 1997 



Wayland, 1997) and glucose-sensing neurons (Schuit et al., 2001 ) have been 



suggested. There may thus be more complex communications in the pan- 
creas for glucose homeostasis, which are left for further study. Finally, we 
also point out that the mathematical model proposed can be generalized to 
describe cellular interactions in other systems, e.g., neural networks consist- 
ing of excitatory and inhibitory couplings. 
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Tables 



Table 1: Coefficients in the transition rate. Here k x = [g' x + jf + jf + 
g x jfj x )(l + g x jf + jf jf + j2 9 x )~ Y with x denoting a, /3, or 5. Parameters 
, jf , and jf are given in Table [2j 



Coefficient 


Value 




1 


< 


k x 




k x j x 


w s 


k x j x 


<p 


AX 
Jl 


w h 


jfjf 


w 8a 


AX 
J2 


K/35 
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Table 2: Parameters in the coefficients of the transition rate. 



X 


a 


P 


6 


9 X 


tanh(G/29) 


tanh(-G/29) 


tanh(-mG/26) 


Ji 


tanh(Ji/46) 


tanh(-Ji/4e) 
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tanh(J 2 /46) 


tanh(J 2 /4G) 
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Figure Legends 

Figure [TJ. 

Schematic diagram of cellular interactions between a-, /?-, and 5-cells. The 
arrow represents enhancement while bars represent inhibition. Here the 
intercellular interactions between two cells are present only when both are 
in active states. 

Figure [2J 

Stationary glucose level Goo, depending on the cellular interaction strengths 
J\ and J2. Note that Goo measures the stationary level relative to the fasting 
glucose level Go in the absence of cellular interactions. 

Figure [3]. 

Time evolution of the glucose level G, depending on the interactions between 
a- and /3-cells. Starting from the high-glucose state [(<r a ) = —1 and (erg) = 1 
at G = ImM (or G = Go + ImM)], the system recovers eventually the 
normal glucose level Goo. The time constants are taken to be r = tq = 1 
for simplicity and the cellular interactions have strengths J\ = 2 mM and 
J 2 = mM. 

Figure 31 

Balance function for glucose regulation, depending on the interactions be- 
tween a- and /3-cells. The strengths of cellular interactions are J\ = 2mM 
and J 2 = OmM. 
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Figure [5]. 

Time evolution of the average states (a a ) and (ap) of a- and /3-cells, starting 
initially from the high-glucose state (a a ) = — 1 and {a/3) = 1 at G = ImM, 
during the recovery to the normal glucose level. The asymmetric interactions 
between a- and /3-cells have the strength J\ = 2 mM whereas the inhibitory 
interactions of 5-cells are absent in (a) J2 = OmM but present in (b) J2 = 
2mM. 

Figure [U 

The average state (cr a ) of a-cells versus (ap) of /3-cells for the asymmetric 
interaction J\ = 2mM, in the absence (J2 = OmM) and presence (J2 = 
2 mM) of the inhibitory interaction of (5-cells. The dotted line along the 
diagonal represents the stationary condition (a a ) = (073}. 

Figure 0. 

Glucose dose-responses in the activity of /3-cells for the inhibitory interac- 
tion J2 = and 2 mM. The asymmetric interactions are taken to have the 
strength J\ = 2mM. 
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